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Porous nanowires (NWs) with tunable thermal conductance are examined as a candidate for 
thermoelectric (TE) devices with high efficiency (ZT). Thermal conductance of porous Si and Ge 
NWs is calculated using the complete phonon dispersion obtained from a modified valence force 
field (MVFF) model. The presence of holes in the wires break the crystal symmetry which leads to 
the reduction in ballistic thermal conductance {(ti). [100] Si and Ge NWs show similar percentage 
reduction in ai for the same amount of porosity. A 4nm x 4nm Si (Ge) NW shows ^ 30% (29%) 
reduction in ai for a hole of radius 0.8nm. The model predicts an anisotropic reduction in ai in 
SiNWs, with [111] showing the maximum reduction followed by [100] and [110] for a similar hole 
radius. The reduction in ai is attributed to phonon localization and anisotropic mode reduction. 



I. INTRODUCTION 

Extreme geometrical confinement and surfaces make 
nanowires (NWs) a promising candidate to obtain higher 
ZT thermoelectric material by suppressing the lattice 
thermal conductivity {k.i). Recent experimental works 
[U [2] reveal that nanowire geometry can greatly enhance 
the ZT of Si (upto 1 at 200K [2 ) from its bulk value of 
0.06 at 300K [3 . An interesting way to further reduce 
the thermal conductivity in bulk Si was shown both ex- 
perimentally [4] and theoretically [5 by creating pores 
in the crystalline material. More recently the fabrication 
of phonon 'nano-mesh' [T showed a 100 fold reduction 
in ni from the bulk Si value of 148 W/m-K to 1.9W/m- 
K. Recent technological developments allow to fabricate 
hollow nanowire arrays using sacrificial templates [6 , hol- 
low spinel wires using the 'Kirkendall effect' [7 , electro- 
chemical anodic dissolution [8 , template based hollow 
wire fabrication [9 , etc. Thus theoretical investigation 
can guide these improved fabrication methods to obtain 
better TE materials and structures. 

The stronger scattering of the phonons from the 
surfaces compared to the electrons, reduces ni drasti- 
cally without affecting the electrical conductivity (G) 
in porous materials due to the much smaller coherent 
phonon wavelengths than electron [5]. Motivated by 
the interesting experimental works and theoretical argu- 
ments, we analyze the ballistic thermal conductance (<j/) 
in hollow Si and Ge nanowires with various channel orien- 
tations ([100], [110] and [111]) to explore the candidates 
which lead to higher reduction in cri and to understand 
the reasons for such reduction. 

The following section (Sec. [ll|) briefly discusses the 
approach for the calculation of the phonon dispersion and 
(Ji in hollow Si and Ge NWs. The effect of the hole radius 
on (Ji a nd t he reasons for the reduction of cji are discussed 
in Sec. Ill, followed by the conclusions (Sec. IV). 
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FIG. 1. (a) Projected hollow Si unitceh with (100) chan- 
nel orientation. Phonon dispersion obtained from the MVFF 
model for (100) (b) hoUow and (c) filled SiNW with cross- 
section size of 3nm(W) x 3nm(H). 



II. THEORY AND APPROACH 



The phonon dispersion in free-standing Si and Ge 
nanowires is calculated using a frozen phonon approach, 
called the modified valence force (MVFF) model [TOHH] 
(Figjlja). In the MVFF method the frequencies of se- 
lected phonon modes are calculated from the forces act- 
ing on atoms produced by finite and periodic displace- 
ments of the atoms in an otherwise perfect crystal at equi- 
librium. The MVFF model has successfully explained the 
bulk phonon dispersions as well as other lattice proper- 
ties successfully in bulk Si, Ge, etc. [TOl [11] and has 
recently been applied to SiNWs to explain their phonon 
spectra and lattice thermal properties [12]. 

From the calculated phonon dispersion, the ballis- 
tic thermal conductance {(Ji) across a semiconductor 
slab/wire can be calculated using the Landauer's formula 
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[13] using Eq.g [HITS], 



(1) 



where M(E), /c^, ^ and e is the number of modes at en- 
ergy E, the Boltzmann's constant, the Planck's constant 
and the electronic charge, respectively. This equation is 
valid when the two contacts across the material slab are 
maintained under a small temperature gradient AT [15] . 

Device details: Rectangular NWs are studied with 
width (W) and height (H) varying from 3nm to 4nm with 
three channel orientations of [100], [110] and [111]. The 
hole radius {vh) in these NWs have been chosen such 
that not more than 25% of the total atoms (Fig. [l] a) 
have been removed from the unitcell to ensure structural 
stability of these wires (no negative phonon dispersion 
is obtained [16]). The inner and outer surface atoms in 
these NWs are allowed to vibrate freely (Fig. [l]a). In 
extremely small SiNWs (W < Inm) significant atomic 
reconstruction can take place which has not been consid- 
ered in this study [17^ . 



III. RESULTS AND DISCUSSION 

Phonon spectra: Hollow nanowires have fewer atoms 
per unitcell compared to the solid nanowires which re- 
sults in a reduced number of phonon sub-bands. Figure 
[l] (b) and (c) compare the phonon spectra of a 3nm x 
3nm [100] SiNW with Vh = 0.8nm and Onm, respectively. 
In the hollow NW a lot of sub-branches appear in the 
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FIG. 2. Ballistic ai in small cross-section rectangular 
nanowires. (a) [100] SiNW, (b) [100] GeNW, (c) [110] SiNW 
and (d) [111] SiNW. The percentage reduction in ai for all 
the wires for ru — 0.8nm is also indicated. 



lower portion of the phonon spectra (Fig. [T] b) indicat- 
ing more surface phonon confinement due to an increased 
surface to volume ratio (SVR). 

Ballistic thermal conductance ai: In all the nanowires 
ai is calculated using Eq.Q. As a general trend, 
ai reduces in hollow nanowires compared to the solid 
nanowires (Fig. |2|. Comparison of ai in Si and Ge NW 
show a similar amount of reduction in thermal conduc- 
tance (Fig. [2] a and b). This indicates a weak material 
dependence of the reduction mechanisms. A 3nm x 3nm 
Si(Ge) NW show a reduction of -37%(38%) for rn = 
0.8nm. This reduction decreases as the wire cross-section 
size increases due to the reducing SVR. The phonons feel 
the surface less in larger cross-section NWs compared to 
the smaller wires. 

The ai reduction is wire orientation dependent as re- 
vealed in Fig. |2]a, c and d. For a 3nm x 3nm SiNW with 
rh = O.Snm, a [100] wire shows a reduction of ~38% 
(Fig. [2] a), a [110] wire shows a reduction of ~27% 
(Fig. 2 c) and a [111] wire shows a reduction of ~40% 
(Fig. 2 d). The ai reduction shows the following order 
[111] ^ [100] > [110] for all the wire cross-section sizes 
considered here. Thus, ai can be tuned by three ways, (i) 
wire cross-section size, (ii) hole radius and (iii) channel 
orientation. 

Reasons for the reduction in ai : The observed trends 
in the reduction of ai can be attributed to two reasons, 

(i) increased phonon localization in hollow nanowires and 

(ii) anisotropic modes and mode reduction in nanowires. 
As the nanowire size reduces the geometrical confinement 
increases which results in increased phonon confinement 
[18 . All the phonon modes in hollow nanowires do not 
propagate well and become localized which reduces ai 
compared to the sohd Si and Ge NWs [T81IT9] . 

(a) Phonon localization: The extent of localization of 
phonon modes can be calculated using a 'participation 
ratio' (PR) [19^. This ratio can be calculated for each 
phonon mode as [19], 

PR-^ = A^S^(Sn,,e[x,^,.]V^^*n,,V^^,n,,)^ (2) 

where N is the total number of atoms in the unitcell, n 
and j represent the number of sub-bands and directional 
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FIG. 3. Average participation ratio (PR) in Si and Ge NWs. 
PR reduces with increasing hole radius indicating increasing 
phonon localization in hollow NWs. 
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FIG. 4. (a) Total modes in hollow SiNWs for three channel 
orientations with hole radius. [110] has maximum modes. 
[Ill] and [100] are very close. (b) Modes distribution in energy 
for three channel orientations. Both plots are for 3nm x 3nm 
SiNWs. 



vectors, respectively and '0i,n,j is the eigen vector asso- 
ciated with atom 'i', sub-band 'n' and direction 'j'. The 
eigen vectors are calculated from the MVFF dynamical 
matrix formulation as described in Ref.[TT]. The partici- 
pation ratio measures the fraction of atoms participating 
in a mode and hence varies between 0(1) for delocal- 
ized states to 0{1/N) for localized states and effectively 
indicates the fraction of atoms participating in a given 
mode. 

The average PR in 3nm x 3nm, [100] Si and Ge NW 
shows a reduction with increasing (Fig. [3]). The 
solid NWs have PR close to 0.7, showing delocalization, 
which gradually decreases indicating phonon localization. 
Thus, both Si and Ge show similar phonon localization 
and hence similar ai reduction (weak material depen- 
dence) (Fig. and b). 



(h) Anisotropic phonon modes: The reduction in cji 
with channel orientation shows anisotropy due to dif- 
ferent propagating modes (M(E)) for each orientation. 
Figure [4] (a) shows that [110] wires have highest modes 
while [100] and [111] wires have fewer modes [14 . [100] 
and [llllwires show similar decrease in total M(E) with 
(Fig- Rla) which explains the similar reduction in ai. 
The energy resolved M(E) for 3nm x 3nm SiNWs with 
Th = 0.4nm (Fig. [4] b) clearly shows higher number of 
modes in [110] orientation thus resulting in smaller reduc- 
tion in cFi. Phonon localization also reflects that a [111] 
SiNW shows similar localization as [100] SiNW whereas, 
[110] SiNW shows less localization (Fig. |3]). This further 
corroborates the anisotropic reduction of cr/ in SiNWs. 

IV. CONCLUSION 

It has been shown that the presence of holes in Si and 
Ge NWs can be used for tuning their thermal conduc- 
tance. Increased phonon confinement, phonon localiza- 
tion due to increased surface-to-volume ratio and mode 
reduction are the reasons for such drastic reduction in ai . 
Thus, variation of nanowire cross-section size, hole radius 
and channel orientation provide attractive ways to tune 
the thermal conductance. [100] and [111] nanowires show 
maximum reduction in gi for Si. Similar trends are also 
expected in GeNWs. This can pave the way to make 
better TE devices using these nanowires. 
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